CMB temperature anisotropies from third order gravitational perturbations 
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' In this paper we present a complete computation of the Cosmic Microwave Background (CMB) 

, anisotropies up to third order from gravitational perturbations accounting for scalar, vector and 

■ tensor perturbations. We then specify our results to the large scale limit, providing the evolution 

' of the gravitational potentials in a flat universe filled with matter and cosmological constant which 

characterizes the Integrated Sachs- Wolfe effect. As a byproduct in the large scale approximation we 
are able to give non-perturbative solutions for the photon geodesic equations. Our results are the 
first step to provide a complete theoretical prediction for cubic non-linearities which are particularly 
relevant for characterizing the level of non-Gaussianity in the CMB through the detection of the 
four-point angular connected correlation function (trispectrum) . For this purpose we also allow for 
O ' generic initial conditions due to primordial non-Gaussianity. 

C/3 , PACS numbers: 98.80.Cq DFPD-A/07/14 
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I. INTRODUCTION 



(N 
> 

00 The three year data set of the Wilkinson Anisotropy Probe (WMAP) on Cosmic Microwave Background (CMB) 

temperature anisotropies has offered a wealth of information about the evolution of the universe of unprecedented 
• accuracy [1] and prepared the way for even more ambitious future missions such as those of the Planck satellite [2] 
[ and CMB polarization observations. Such precise measurements have been accompanied by increasing theoretical 
C , efforts in predicting the signatures on the CMB from various cosmological scenarios. In particular a lot of attention 
in the last years has been dedicated to the statistical properties of the CMB beyond the power spectrum, in search for 
possible non-Gaussian signatures [3]. Non-Gaussianity of the CMB owes its importance to the possibility of unveiling 
crucial aspects of the physics of both the early and the late universe, which would be unreachable using only the 
rS CMB power spectra information. Different mechanisms for the generation of the cosmological perturbations predict 
different amplitudes and shapes of primordial non-Gaussianity, thus a positive detection of NG or an upper limit on 
its amplitude is powerful in discriminating among the various competing scenarios which would be indistinguishable 
otherwise [3] . An illuminating example is given by the standard models of inflation which typically predict a very 
low content of non-Gaussianity [4, 5] and as such they could be completely ruled out if a positive detection will be 
achieved. CMB non-Gaussianity can also have a non-primordial origin due to secondary anisotropies which arise when 
the CMB photons leave the last scattering surface and cross the large-scale structure of the Universe. These include 
both secondary scatterings (such as the thermal and kinetic Sunyaev-Zel' dovich effects, produced by the thermal 
and bulk motions of electrons in clusters, and the Ostriker-Vishniac effect, due to bulk motions modulated by linear 
density perturbations) and gravitational secondaries such as the gravitational lensing and the non-linear Integrated 



'Electronic address: damlcoasissa. it 
tElectronic address: nicola.bartoloOpd. infn. it 
Electronic address: sabino.matarreseapd.infn.it 
^Electronic address: antonio.riottoapd.infn.it 



2 



Sachs- Wolfe effect or, on smaller scales, the Rees-Sciama effect due to the non-linear evolution of the potentials. Such 
non-linearities can give important and new information about the dark matter and the dark energy content of the 
universe [6-13]. Statistics like the bispectrum (the Fourier transform of the three-point correlation function) and the 
trispectrum (four-point correlation) can then be used to assess the level of primordial non-Gaussianity on various 
cosmological scales and to discriminate it from the ones induced by the secondary anisotropics and systematic effects. 

Up to now most of the attention has focused on the three-point statistics both for the primordial non-Gaussianity 
and for the analysis of the angular bispectrum of the CMB temperature and polarization anisotropies [3, 10, 14-17]. 
On the other hand, the four-point correlation function could display interesting features. From the observational point 
of view it has been argued that it could be even more sensible to primordial non-Gaussianity than the bispectrum 
for very small angular scales in the next CMB experiments [18, 19]. ^From the theoretical point of view the four- 
point connected correlation function for the large scale CMB anisotropies has been computed in Ref. [20] giving 
the theoretical predictions for the quadratic and cubic non-linearities (which include generic non-Gaussian initial 
conditions), while Refs. [21, 22] have computed the trispectrum from single-field and uncoupled multiple fields in 
slow-roll inflation showing that it is of the order of the slow-roll parameters. Interestingly in some cases, such as in 
some configurations of the curvaton scenarios [23, 24] , the main source of a non-Gaussian signal can come from the 
four-point correlations functions, the three-point correlation function being suppressed (see also the phenomenological 
model discussed in Ref. [1] and Ref. [25]). 

However most of these computations deal with the trispectrum of the uniform density curvature perturbation ( on 
large scales within a given inflationary model, but this is not the physical quantity that is observed. The same is 
true of course for the bispectrum of the curvature perturbation. But, if for the case of the curvature perturbation 
bispectrum one is sure that the information about the primordial non-Gaussianity in the flnal observable quantity 
(the CMB anisotropies) is just the one obtained by evolving linearly the curvature bispectrum, the same is not true 
any more for the trispectrum. Let us consider the standard way to characterize the level of non-Gaussianity in the 
gravitational potential by expanding it as 

$ = *L + /NL($i-(*i))+flNL$i, (1.1) 

where $l is the linear first-order Gaussian part, /nl and ^nl arc the parameters which measure the quadratic and 
cubic non-linearities. In fact the non-linearity parameters might have in general a non trivial scale dependence. At 
linear order the curvature perturbation is ("'^^ = 5$l/3 during the matter dominated epoch, and, for example, the 
Sachs-Wolfe effect tells us that AT/T = — $i,/3. The standard approach followed in the literature is to extend both 
of these two relations also at higher-order, so that for example for quadratic non-Gaussianity one writes C = 5<I>/3 
using the expansion (1.1) up to second-order. Evolving linearly the perturbations in such a way just accounts for 
the primordial content of non-Gaussianity and is justified only when one assumes that the primordial level of non- 
Gaussianity is much larger of the second-order corrections (in the form of first-order squared perturbations) which 
arise both in the relation between the curvature perturbation and the gravitational potential and when computing 
the CMB anisotropies. One is guaranteed that such primordial contribution in the final CMB anisotropies will appear 
as computed in this way because it is already an intrinsically second-order non-linearity and all the transfer functions 
will therefore be the same as at linear order. For the CMB bispectrum this is also all what is necessary to account 
for the primordial non-Gaussianity. For the trispectrum the situation is different, however. Evolving linearly the 
primordial cubic non-Gaussianity in Eq. (1.1) would just lead to (^^^ /6 = 5(7nl/3$|, (splitting ( = (^^^ + (^(^'/2 + ... 
into a first- and higher order parts). Since the curvature perturbation remains constant on superhorizon scales (for 
adiabatic perturbations) it is a useful quantity to keep track of the primordial non-Gaussianity. Therefore, suppose 
to parametrize the primordial non-Gaussianity as 

C = C^') + Kl - 1) (C^^^)' + (6nl - 1) (C^^^)' , (1.2) 

where the two non-linearity parameters depend on the physics of a given scenario for the generation of the perturbations 
(for example, for standard single-field models of infiation cnl = 1 and 6nl = 1, plus tiny corrections proportional to 
the slow-roll parameters). Thus the relation C^^V^ = Sanl/S^I, is equivalent to 

25 

fifNL = y(&NL - 1) . (1.3) 

However, this does not catch at all the whole information about the primordial non-Gaussianity that is contained 
in the CMB anisotropies up to third order. This is due to two reasons: first, as a source for the evolution of the 
gravitational potentials at third-order now there are also the gravitational potentials at second-order which inevitably 
contain primordial non-linearities proportional to the /nl (or cnl) parameter; and second, when computing the CMB 
anisotropies at third-order, one still gets additional terms proportional to /nl^l- This point is clear in the results 
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we present in Sec. Ill and IV. See for example the first line of Eq. (3.5), or Eqs. (4.2) and (4.7). Another example 
can be found in Rcf. [20] , where a fully non- linear expression for the Sachs- Wolfe effect has been obtained for generic 
non-Ganssian initial conditions. The corresponding parameter (/nl for CMB anisotropics has been computed (see 
Eq. (67) of that Reference) by accounting for all the possible dependence on the primordial non- Gaussian! ty finding 
an expression of the type 

25 

5nl = yC'NL - 1) + A(ki,k2,k3;aNL), (1-4) 

which, contrary to Eq. (1.3), depends also on the quadratic non-linearity parameter (here kj are the perturbation 

wavenumbers in Fourier space). A similar expression to Eq. (1.4) is given for the gravitational potentials at third-order 
in Eq. (4.23). These examples show that even in the case of a "local" model for quadratic non-Gaussianity (i.e. a 
constant /nl parameter), at higher orders such non-linearities are modulated by first-order perturbations generating 
a scale-dependent non-Gaussianity. 

The best limits up to date on the /nl parameter come from measurements of the CMB bispectrum on the WMAP 
data giving —36 < /nl < 100 [1, 14, 26], while there is at present no real bound on ^nl- Given the increasing 
precision of future mesurements of CMB anisotropies, it is clear that it is of fundamental importance to provide 
accurate predictions for all the cubic non-linearities that enter in the evolution of the cosmological perturbations. 
This is not only mandatory to be able to evaluate the trispectrum of CMB anisotropies, but becomes crucial if 
one wants to account for the precise dependence of the trispectrum on the primordial non-Gaussianity. Spurred 
by these considerations, in this paper we will focus on the CMB anisotropies up to third-order from gravitational 
perturbations due to the redshift the photons suffer when they travel from the last scattering surface to the observer. 
The computation will be performed following Refs. [27-29] by perturbing at the desired order the photon geodesic 
accounting for scalar, vector and tensor perturbations. However we will also propose an alternative method to solve 
the photon geodesic equation (for scalar perturbations), which consists in a non-pcrturbativc computation. In the 
last part of the paper we focus on the large-scale approximation and on scalar perturbations for a Universe filled 
with non-relativistic matter and a non-vanishing cosmological constant. In fact one of the goals is to characterize 
the evolution of the gravitational potentials on large scales which is responsible for the (late) Integrated Sachs- Wolfe 
effect, thus completing the analysis for the large scale CMB anisotropies at third-order already started in Ref. [20]. Of 
course it is out of the goals of this paper to deal with the anisotropies generated from the dynamics at recombination 
and due to the scattering terms ^. 



II. TEMPERATURE ANISOTROPIES 



We arc interested in the pattern of fluctuations of the CMB temperature as measured by an observer in a perturbed 
flat Friedmann-Robertson- Walker spacetime. 
The line element can be written as 



„(o) 



.„(l)+„(2) 



.„(3) 



dxi'dx" . 



(2.1) 



where 0(77) is the scale factor in conformal time drj = dt/a{t), g^J is the background Minkowski metric, and g'j^J is 
the r-th order metric perturbation. 

Photons travel along null geodesies x^{X), where A is an affine parameter in the conformal metric g^^. The photon 
path connects the point of observation, with coordinates Xq = {i]o,xo), to the hypersurface of emission, deflned as 
the spacelikc hypersurface of constant conformal time 77^ . The actual last scattering surface is the intersection of the 
observer's past light cone with this hypersurface. We assume that on the constant rjg hypersurface every point p 
emits thermal radiation, characterized by a temperature r(p, d) which depends on the point of emission p and on the 
direction of emission, described by the vector d normalized to unity in the conformal background metric. The different 
photon paths are specified by the direction from which they arrive at O, described by a vector e, normalized to unity 
in the conformal background metric; this unit vector can be thought as the direction toward which the observer is 
pointing an antenna. The initial conditions Xq, e determine the point p and direction d of emission. 



A complete computation of the CMB anisotropies due to the non-linear dynamics taking place at the last scattering epoch has been 
performed in Refs. [30, 31], where the full system of the Boltzmann equations up to second-order for photons, baryons and cold dark 
matter have been presented togheter with analytical solutions in the tight coupling approximation. 
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During their travel from the last scattering surface to the observer the CMB photons suffer a redshift determined by 
the ratio of the emitted frequency los and the observed one loq- For a blackbody spectrum, the ratio uj/T \s constant 
along the photon path, and the temperature measured by an observer is given by 

ro(xo,e) = ^^T£(p,d). (2.2) 

The expression for the frequency is 

w = , (2.3) 

where u'^ is the four- velocity of the observer or emitter, normalized to a^QuvU^u" = —1, and fc'^(A) = ^ is the 
photon wavevector, tangent to the null geodesic x'^{X). 

Given the initial conditions XQ,e,ujo, we need to compute the quantities p,d,ti;£, propagating the photons back 
from the observation to the emission surface. These will depend on the photon path and the associated wavevector, 
which we expand in series of the metric perturbations gj[J and their derivatives: 

a;" (A) =a; (A) + aj^^)" (A) + a;^^)'' (A) + x^^^i" (A) + . . . (2.4) 
fc''(A) =fc(°)''(A) + fc(i)''(A) + fc(2)''(A) + + . . . (2.5) 

We find it useful to write the perturbed (conformal) metric as 

gf.^dx^dx" = -e^*d?7^ + 2u>idr]dx^ + {e'^'^Sij + Xy ) dx^dx^ , (2.6) 

which is valid at any order in perturbation theory. The two gravitational potentials $ and "l/ correspond to scalar 
metric perturbations, uJi includes a scalar and a vector perturbation, and Xij contains another scalar, as well as vector 
and tensor perturbations. To make contact with the usual perturbative calculations one expands every quantity as 
$ = $(1) + _|_ i$(3) _|_ _ _ g^j^^ analogously for the others. Notice that, unlike u'f '' and Xt'j\ the quantities ^^""^ 

and ^^^^ are not the usual r-th order scalars (^^''^ V^'^' appearing in Refs. [3, 32]. However, it is immediate to find the 
relation between these quantities. Up to third order $(1) = ^^^^ = 'ip^'^\ = (^(^^ - 20^^^, ^^^^ = i/'^^^ +2?/'^i), 
$(3) = ,/)(3) - 60(1)0(2) _^ g^3^^^ ^(3) ^ ^^(3) ^ Q^W^C^) + 8^^3^^. The form of the metric (2.6) greatly helps in 

the intermediate computations, however some results will be expressed in the variables 0''") and ^^"^^ because they 
appear in a more compact form. Notice that in this section we will not choose a particular gauge so all the following 
expressions for the CMB anisotropics are valid in any gauge. 
The four- velocity will be expanded as 

u^^ = l(s^ + 1,(1)^' + ^(2)/^ + L(3)/^ + . . . V (2.7) 
a \ 2 6 / 

The zero component of the velocity is fixed from the normalization condition; we find 

^["d = - -^^'^ ; (2.8) 



+ 3(0W) +2u;Pvl^+vl,^v^'^; (2.9) 
v1,^ = - <^^') + Sco^vl^ + 3a;f + 90(^)0(1) + 3 (xg^ - 2V>(i)^,,) v^^vl^ 

- Ucl^^M'^vl) - 15 (0^'^)' - 30(i)z;^i)^;f ) + Sv^^vl'^ . (2.10) 
In order to compute the observed temperature up to third order in perturbations, we expand the frequency as 

w = (1 + + + Co^3) ^ _ ^ (2.11) 

and the temperature at emission as 

T£(p,d)=Tf'(l + T(p,d)) , (2.12) 

where r ~ r^^^ + t^^-* + r^^) + . . . is the intrinsic temperature fluctuation at emission. A calculation of this quantity 
up to third-order is beyond the goal of this paper, since we are interested in the additional effects of gravity along 
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the photon path. However on large scales (bigger than the horizon at recombination) it has been computed in a 
non-pertubative way in Ref. [20]. Its fully non- linear expression is very simple: t = e~^*/'^ — 1. Moreover, on 
smaller scales, a complete treatment of the CMB anisotropics including the acoustic oscillations on the surface of last 
scattering up to second-order have been studied analytically in Ref. [30, 31], and a full numerical evaluation can be 
performed Ref. [33] using the set of Boltzmann equations provided in Ref. [30] . 

Now, we have to take into account that we need to expand the point p and direction d at emission as p = 
p(°) + p(^) + p^^^ + . . . and d = d^^^ + d^^' + d^^^ + . . . ; then we Taylor expand t(p, d) around the background values, 
obtaining 



^(P,d) 



+ d} 



r(3) 



^2^(1) 



(2.13) 



1 , J 



(2) dd^ 



where it is understood that all quantities are evaluated at p*-"'*, d*'"-*. We find p^^^^ = x^^^^ — k^^^^x^^^^, p^^^' = 

^(2)4_ l(0)*^(2)0 , r.(0)zjL(l)0^(l)0_ ^,(1)1^(1)0 _ fc'°''+fc(i)' fc'"'' ^(2)» _ fcW'+fc'D'+fet^)' fcW'+fcd)' t^„4.;„„ 

X fi.x-h«,«,x n X ,a - |^(o)i_^j,(i)i| I^MTj' " - |fe(o)i+fe(i)i+fe(2)i| |fe(o)i+fc(i),|- i^owt,e 

that therefore the quantities d^*^ are just the difference between the photon normalized wavevector up to the i-th 
order and that at the {i — l)-th order. 

Performing all the expansions in Eq.(2.2) and dividing by the background temperature ^T^\ we find 



+ T 



(1) . 



(2.14) 



(2.15) 



' "(^^ 2^(1)^(1) dx^ 
avd) , ar(i) „ arfi) 



, (2) ^ 5^^'^ ,n dT^^'^ -(2) 



dx^ ■ "w^""^^ ; (4'^-4'0 
+ (r(i) - 4^)) - 4^) + (4^V - 4^^4^0 • 

The expansion of the frequency using Eq.(2.3) yields: 

4i) = fcO + <^(i)+4i)e'+^We^ 



(2.16) 



(2.17) 



4'^ =kU + l^''' + \^?'e^ + lvl,v^' - ^(<A«)^ + - 4^^fcb - 



aA 



+ 



(2.18) 



6 



(2.19) 



,(1) 



rffc 



(1) 



"3 "(1) ^ "(1) "3 



^-^'''^^3)+pU^,{^-2^(^^S,J) 



dX 



+ 
+ 



In these expressions, AA^^^ and AA^^^ + AA^^^ are the differences in affine parameter between the points where the 
first order (respectively, the second order) and the background geodesies intersect the r] = rjs hypersurface. They are 
given 

The next step is to obtain the null geodesies up to third order, using the formalism set up in [28, 29]. This will 
allow us to compute at the desired order the photon path a;^(A) and the associated wavevector A;'' (A) appearing in the 
previous expressions for the perturbations of the photon frequency. The a-th order geodesic equation can be recasted 
as the forced Jacobi equation [28, 29] 



rfA2 ' ""'^ 

In a flat background, the solutions are 



+ 2r(°)''fc(°)"fc(")'^ + 5<^r(°)''fc('')"A;(°)'^x('')'^ = . 



=(A- Ao)fc(«)''(Ao) + / (A- A)/(")''(A)rfA, 

^ Xn 



The forcing terms, in a flat background, are given by 

J (1) — al3 ' 



- _ or(l)^l,(0)ai,(l)/3 _ T.W<ya r(l)'^l,(0)ai.(0)/3 _ p(2)M i,(0)a , (0)/3 



/(2) 
/(3) 



1^^'*^(2)'«(0) + ^'^(l)'^(O) + ^-^^ a/3 ^(1) J 

'^(O)'^(O) I +"<^^a0 ^(2) + O'a-L a/3 + 2 " ^(1)^(1 



(2.20) 

(2.21) 
(2.22) 

(2.23) 
(2.24) 

(2.25) 
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The ChristofFcl symbols are calculated in the Appendix, where one can find their expressions for the conformal metric 
(2.6). Notice that we have recastcd them in a compact form without specifying the order of the perturbations, even 
tough they can be used just up to third-order (at least for the terms coming from the vector and tensor perturbations). 
The method that uses Eqs. (2.23) is an iterative method: once one knows the geodesic equations oir — th oder, he or 
she is able to determine those at the next-order through Eqs. (2.21). Thus, in order to have all the means to compute 
the photon geodesies up to third order in the following, we report the expressions that we have computed for the 
wavevectors fc^'*) up to second-order : 



fcO^j = -2$(i) - ^e' + / d~XA^^^ (A) = -2$(i) - ^e' + /i(A) , 

J Ao 



(2.26) 



(2.27) 



for the first-order wavevectors, and 



+ 



(2.28) 



+ :^X(2)ie^+a;^c*^X(i),e^-X(i),- 



+ / dX 



(2.29) 



at second-order, where 



(2.30) 



As a non-trivial test for the correctness of the above expressions, we have explicitly verified that the first and second 
order wavevectors do satisfy the null vector condition k^kp, = 0. When expanded perturbatively, this reads at linear 
order 



and at second order 



fc(2) + e.fc^2) + - + 2$(i)fc(i) - ^4) + '^P^e'^) - 2*(i)e.fc^i) + xg^e'fcf^^ 

Both equations are satisfied when one substitutes the expressions (2.26)-(2.29). 



(2.31) 



(2.32) 



III. CMB ANISOTROPIES ON LARGE SCALES FROM SCALAR PERTURBATIONS 



In this Section, we focus on scalar perturbations and we give the expressions for the CMB anisotropics on large 
scales. Two main results are achieved. First we are able to give, using a non-perturbative method, an integral 
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solution for the photon geodesic equations which holds at any order in perturbation theory. Second, we solve for the 
first time the evolution of the third-order gravitational potentials in a universe filled with non-relativistic matter and 
a cosmological constant. The presence of a non- vanishing cosmological constant makes the gravitational potentials 
vary at late times producing on large scales an Integrated-Sachs- Wolfe effect which is the counterpart of what happens 
at linear order. Of course, due to the non linear evolution of the perturbations, other integrated terms will appear 
in the form of cubic corrections. The ISW constitutes the main effect for large scale CMB anisotropies together with 
the Sachs- Wolfe effect originating at the last scattering surface. In fact both of them keep memory of the initial 
primordial non-Gaussianity of the cosmological perturbations (see details in Sec. IV B), and as such our expressions 
can be of particular relevance when trying to pin down the primordial non-Gaussian content from the higher-order 
statistics of the GMB anisotropies as the bispectrum or the trispectnim. Let us sec in details these two results. 

In the limit of large scales, it is convenient to adopt a non-perturbative formalism to get the temperature anisotropies 
and the Einstein equations. We take the metric to be ds^ = a^{ri)[—e^'^ drj^ + e^'^'^ Sijdx^dx^] and for simplicity we 
take the points of emission and observation to be comoving (in such a way we lose the Doppler effects, however they 
are important on scales smaller than those we are interested here). In writing this metric we have neglected vector 
and tensor perturbation modes. For the vector perturbations the reason is that we are interested in long-wavelength 
perturbations, while vector modes will contain gradient terms being produced as non-linear combination of scalar 
modes and thus they will be more important on small scales (linear vector modes are not generated in standard 
mechanism for comsological perturbations, as inflation). For example the results of Ref. [34] show clearly this for 
second-order perturbations. In order to study the CMB anisotropies from scalar perturbations in the large scale limit, 
the tensor contribution can be negleted, since on large scales it has been proven to remain constant and to give a 
negligible effect being of the order of (powers of) the slow-roll parameters during inflation [5, 35]. 

We start from eq. (2.2) and we write the frequency as 

CO = -g^^u'^k" = e'^^u°k° = e*fc° , (3.1) 
where = e~* is given by the normalization condition. The observed temperature is given by 

ro(x,e) = ^e*--*^§T£(p,d). (3.2) 

We make another simplification: namely, we reabsorb the terms computed at the point of observation into a 
redefinition of Tq. This amounts to lose the monopole term, which however is unobservable. 

At this point, we can perform the double expansion in perturbation orders and around the background geodesic 
(i.e., the line of sight). We find 

^=r«-4^)-fcr; (3.3) 



■d\. 



dd' 



(2) 



dk 



(1)0 



d\ 



+ (fcW0)2_^(l)$(l)_^(l)fc(l)0 + $(l)fc 



(1)0 



(l)^2 
S > 



(3.4) 



1.(1)0. 



5(3) T 



.(3) - 43)0 _ 1^(3) ^ 1^(1)^(2) _ 1 (^(1)^)3 ^ _ ^^(,^„^3 _ ^^^^^^^^ 



+ 



42)0 _ (^(l)0)2j (^$(1)_^(1))+V) 

a.^(2) _ i5.$(2) _ ^(i)a.$(i) _ $(i)a,r(i) + (dM^^ - dM'^) + 



■P(i) 



(1) 



dd' 



dx^dd^ 



dX 2 dA2 



(1)0 



+ (2fc«° + $«-rW)^^ 



,dk 



(1)0 



-aa(2) ^ . 

dX 



(3.5) 
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In order to compute the photon wavevectors and the null geodesies, we will proceed in a different way than the 
usual standard perturbative method introduced in [28, 29] which makes use of the Eqs. (2.23)-(2.25). Let us instead 
start from a fully non-linear geodesic equation (time component) which can be written as 



(3.6) 



where we have used the null vector condition g^i^k^k'^ = 0. We can write the formal solution (formal since the 
potentials depend on the true geodesic) as 



fcO(A) = e-2*(^) 



-2*(x(A)) 



($'(a;(A)) + *'(a;(A))) d~X 



(3.7) 



The advantage of Eq. (3.7) is that it can be straightforwardly expanded to find the perturbations to the background 
wavevector (and, by integration, to the background geodesic), up to any desired order. To make contact with the 
standard perturbative results, in addition to the expansion in perturbative orders we need to perform a Taylor 
expansion of the true geodesic around the background geodesic for the potentials, in order to express all the quantities 
on the background geodesic. For example, if we have the non-linear quantity a{x) computed on the true geodesic, its 
expansion up to second order will be a{x) = a^^\x^^^) + di^a^^\x^^^) x^^^^ + a^'^'{x^^^), and similarly at higher-orders. 

We find 



A;(i)0(A) = -2$(i)+7i(A), 



(3.8) 



fc(2)0(A) = -$(2) + 2 - 2a;(i)''9^$(i) + 72(A) - 2$(i)7i(A) + h{Xf , 



(3.9) 



where 



A;(3)0(A) = - _^(1)A'5^^(2) _ 2x(2)''a^$(l) + - ($(^))' 

3 3 ' 

+ (2$(i) - 7i(A)) [$(2) - 2 + 2x(i)^a^$(i) 

+ 73(A) + 27i(A)72(A) + h{Xf - 2$(i) (72(A) + 7i(A)2) , 



/i(A) 
72(A) 
/3(A) 



Xo 
A 



dX 



1 A(2)'(A) + x(i)^a^A(i)'(A) - 2$(i)A(i)'(A) 



dX 



^' 1 



- 2$(i) i^u^^y + x^'^'d^A^^yj + A^^y (2$(i)' - $(2) _ 2x(i)^a^$(i) 

j^{n) =$(n) ^ ^(n) _ 



(3.10) 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



The above expressions arc completely known once we find the geodesic x^{X) (it is sufficient to compute them up to 
second-order, as it is clear by looking at Eq.(3.13)). For the time component the integration of fc*^°)(A) yields 



c(i)0(A) = dX [-2$(i) + (a - a) 



(3.15) 



a;(2)°(A) = 



dX -$(2) + 2 ' - 2x(i)''a^$(i) - 2cl>(i)7i(A) + 7i(A)2 

+ (a-a) QA(2)'(A) + a;«''a^A«'(A)-2$WA«'(A) 



(3.16) 
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Notice that from Eq.(3.6) to Eq.(3.13) wc have never used the large scale approximation, so these expressions are 
general and we recover exactly the corresponding expressions for the second-order quantities already found in Refs. [27, 
29]. On the other hand to compute the spatial components x^{X) we do make use of the large scale approximation in 
the following way. We start from the non-perturbative spatial geodesic equation 



dX dX ^ ' 



(3.17) 



To find the solution of this equation we first decompose fc* = fc||e* + k\ where fcjj — e^W and k\ — ((5* — e^ej)k^ are 
the parallel and transverse parts of the wavevector with respect to the background geodesic. Then, in the geodesic 
equation we approximate S.^jk^k'^ = fcj^ + Sju^^k^ ~ fc^, since the second term will be negligible on large angular 
scales. Then, projecting Eq.(3.17), we can split it in the two equations 



dX 
dk 



dX 
d^a 



which can be solved. The solution for kn is analogous to the one for fc" 



(3.18) 
(3.19) 



fell 



„2* 



-1 



/ dXe'^^e'di (* + $) 



from which the solution for kj_ is easily obtained 



ki = -e'* / dAe'*((5*^' - e'e')dj (* + $) -1 + / dXe'di (* + 

J Xo Xo 



(3.20) 



(3.21) 



since as initial condition we have k\{Xo) = 0. Perturbing Eqs. (3.20) and (3.21) up to second-order we find 



^(1 



(1) 



.2^-(i)e* - / dXd'A'-^'^ 



(3.22) 



+ f dX 

Xo 



.^5^(2) _ - 2*(i)aM(i) - 25^(1) / dX'e^djAC"^ 



(3.23) 



Notice that in fact the above expressions for fc' are exact, as we have checked explicitly by using the pertubative method 
of Eqs. (2.23)-(2.25). The reason is that up to this order the large-scale approximation we used also corresponds to 
neglect higher-order terms in the perturbations. Moreover, as a further consistency check, we have explicitly verified 
that the expressions (3.8)-(3.10) and (3.22)-(3.23) are the same one gets with the perturbative method employed in 
Sec II. 

The integrals (3.11)-(3.13) along the line of sight are at the origin of the late (and early) Integrated Sachs-Wolfe 
eff'ect, since they express the redshift the photons suffer from the last scattering surface to the observer while traveling 
through a time varying gravitational potential. In particular the first integrand terms, ^^^^ and A^^^^ in each of the 
integrals at second and third order, /2(A) and /3(A), are just the straightforward extension of the well-known first 
order ISW effect. However one must consider all the additional contributions which are integrated terms and which 
depend on the time variation of the gravitational potentials. So terms of this type will appear also from the first two 
lines of Eq.(3.10). Notice in particular that terms containing spatial gradients like x^^^^djA'-'"^ (and similar) coming 
from the combination x^'^'^^dfj^A^"^'^ = x^'^'^^ djA^"^'' + x^'^'''^doA^™^ cannot be neglected a priori even in the large scale 
approximation because they are integrated along the line of sight. ^ 



^ In fact one can check that for a term Uke, for example, x^^^^ djA^'^'^ only the parts of such gradients transverse to the background 
geodesic can be neglected on large scales, while the longitudinal part combines with x^^^'^doA^"^^ to give a term which does not contain 
spatial gradients. 
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Finally plugging the expressions (3.8)-(3.10) for k^^^ up to third-order into Eq.(3.5) the expression for the third 
order temperature fluctuations on large scales is 



= r(3) + + l^'M' + I 



+ x 



he — he 

(1)0 



(2) + ^(i)$(i) + ^ 1^(1)^ ^ 2x(i)°4' 



1 



T^^s ^2 e 



9rW 



+ heAf + he (-24^)' + 4^0 " r^'"' K'" ' 
- (2he - 3$^'^ - r(i)) - 2<^f) ] + 2AA(2)$«' - AA^^)^)' , 



(3.24) 



fc'-^'o.T^^^*'. Wc recognize various contributions to 



where ^^""^ is given by Eq. (3.14) and wc recall that AA*^^) — — x^^'" - 
the temperature anistropies in Eq. (3.24). The first line includes the Sachs- Wolfe effect, which combines the intrinsic 
temperature fluctuations t with the photon redshift on the last scattering surface due to the gravitational potential 
perturbations The terms in the third line correspond to the Integrated Sachs- Wolfe effect, while in the fourth 
line there are contributions due to a lensing effect at the last scattering surface, being dependent on the vector d 
which specifies the direction of emission of the photon. Finally all the remaining terms are due to a possible time 
dependence of the gravitational potentials at the last scattering epoch, due to the fact that by that time the universe 
is still not completely matter dominated (similarly to the early integrated Sachs- Wolfe effect). At second-order, in 
Eq. (3.4), the lensing like contributions appear as well. Notice that in order to make a comparison with the bispectrum 
computed using the different technique of Ref. [36] , one has to take into account also these lensing terms along with 
the Sachs- Wolfe effect. 



IV. THIRD-ORDER SCALAR PERTURBATIONS OF A FLAT ACDM UNIVERSE 
A. Evolution of the gravitational potentials on large scales 

We now consider a spatially flat Universe filled with a cosmological constant A and a non-relativistic perfect 
pressureless fluid and we solve the third order Einstein equations for the scalar perturbations, in the large scale limit. 
The energy momentum tensor of the matter component reads Tj^ = pu^^u^, with energy density p and four-velocity 
w'*. The derivation of the relevant equations is sketched in the Appendix. 

The evolution equation for the third order gravitational potential V''^'^^ reads 

Vfs) + 3WV(3) + «'AV'(3) = S{71, x) , (4.1) 

where the source term is given by 

5(r?,x) =^^±-V-2{ [{n' + a'K)E + nE'] [2,V-''^^^' ((^oaVo9j<^o) - (^o(5Vo)(9i<^o)] 

+ [{H' + a^A) F + TLF'] [iV-'^d^d' (d'aodjipo + d'ipodjao) - 2{d^ao){di^o)] 3) 

- ig' + 2ng) [m{g'f + a^AgH^f - 1)] ^o^o} 

- 6(5') Vo + 3 {g'A'ifl + g'B'ifoao) + gradients , 

where stands for the inverse of the Laplacian operator. Here g{r]) is the linear growth function, defined by 
(p{'x.,ri) = g(?7)(/?o(x), where 9?o(x) is the peculiar gravitational potential linearly extrapolated to the present time, 

and f{r]) = 1 + The growth-suppression factor is given by g{r]) = D+{r])/a{ri) where D+{ri) is the linear 
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growing-mode of density fluctuations in the Newtonian limit. The exact form of g can be found in Refs. [37-39]. A 
very good approximation for g as a function of redshift z is given in Refs. [37, 38] 

gocQm [flT - ^^A + (1 + 0^/2) (1 + fiA/70)] ' , (4.3) 

with n„, = + zf/E^{z), 17a = n^K/E^iz), E{z) = (l + z)H{z)/no = [l7o™(l + zf + floA]'^^ and Oo™, 

l^oA = 1 ~ l^Om, the present-day density parameters of non-relativistic matter and cosmological constant, respectively. 
We will normalize the growth-suppression factor so that g{z = 0) = 1. The function /(ry) can be written as a function 
of Qrn as /{n^) « n^iz)^/"^ [37, 38]. In the A = case 5 = 1 and /(r/) = 1. 

In Eq.(4.2) A{r]), B{r]) enter in the solutions for the second order gravitational potentials on large scales [13]: 

*(2)(x,r,) =A{r,)^l{^)+B{r,)ao{^) (4.4) 

$(2)(x,r,) =AM{^) + C(r?)ao(x) (4.5) 

and their explicit expressions are 

ao(x) =V-2 [aVoa.</^o - 3V-29,9J(9Vo d^^o)] , (4.6) 

A{r,) =Bi{tj) + 2g^{7j) - 2gmg{v) - y (ajVL - 1) 9mg{v) , (4-7) 

B{rj) =B2{v) - IgmQiv) , (4.8) 

C{v) =B{v) + ^gmgiv) (P^ + I) , (4.9) 



where the fimctions -61(77) and B2{f]) arc given in Ref. [13] while g{f]rn) = 5m is the value of the growth suppression 
factor during matter domination, when the cosmological constant was still negligible. A good approximation is 

Notice that in A{r]) appears the parameter a^L specifying the level of quadratic primordial non-Gaussianity which 
depends on the particular scenario for the generation of cosmological perturbations, as we will discuss later in more 
detail. For example, for standard single-field models of slow-roll inflation cnl = l-l-0(e, rj) [3-5, 40], where e and r] are 
the standard slow-roll parameters [41]. The reason it appears in Eq. (4.2) is due to the fact that at third-order some 
of the source terms contain products of linear and second-order gravitational potentials, as detailed in the Appendix. 

We have also introduced the functions 



/\2 



E{r,) =2A'g' + {bH^ - a^A) Ag - 2ng'g^ + 2n{A'g + g'A) - 2g{g' 

- 2g{g' + Hgf + 2^ {g' + Hgf , ^^'^"^ 

Fin) ^^B'g' + \ - a' A) Cg + \n{Cg' + B' g) . (4.11) 

Notice that when A = then Bi{'q) and 52(7?) 0,B(??) -4/3 and A{ri) -^2- (10/3)aNL, and therefore Eijf) 
simplifies to E{ri) = -(50/3)aNLW^ while F{j]) = (5/2)^2. Thus, looking at Eq.(4.2), it is simple to see that all the 
terms in the last two lines are specific only to the case of a non-vanishing cosmological constant. 
The solution of Eq.(4.1) is given by the Green's formula 

i^^'Hv) = jj^^^ + i'+iv) £ dr?'|^5(V) + ^_(r;) £ dv|±M5(r;') , (4.12) 

where 

V^+(^) = g{^) , V-(ry) = ^ , W{rj) = ^ (/o + l^^m) , (4.13) 

are respectively the growing and decaying mode solutions and the Wronskian of the homogeneous equation. Here the 
suffix '0' stands for the value of the corresponding quantities at the present time, while = 4'''^\ilm) represents 
the initial condition taken deep in the matter dominated era on super-horizon scales, r]m being the epoch when full 
matter domination starts. It is such an initial value that must be properly determined in order to account for the 
primordial cubic non-Gaussianity in the cosmological perturbations. 
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The solution for cj)^^^ is then obtained from the relation between tp^^^ and ^'^^ obtained in the Appendix (Eq. B7) 

^(3) = ^(3) + i2g(r^)A{r^)^l + 6g{v) {B{rj) + C(r;)) ao^o - [E{f^)fio + F{t^)vo] , (4.14) 

where we have defined the following kernels: 

/xo(x) =3V-%d^ (<^o 5Vo dj^o) - (<^o ^fe^o) , (4.15) 
ivo(x) =3V-%d^ {d'ao dj^o + d'^odjao) - 2y-\d''ao dk^o) ■ (4.16) 



B. Initial conditions from primordial non-Gaussianity 

We now discuss the key issue of the initial conditions, which are conveniently fixed at the time when the relevant 
modes of the perturbations are well outside the Hubble radius. In order to follow the super-horizon evolution of the 
density perturbations produced during inflation, we use the curvature perturbations on uniform- density hypersurfaces 

which will be expanded as C = C''"^^ + sC*'^^ + gC^^'' + • ■ • ; where at linear order C*^-"-) = —tp^^'> — HS^^^p/p was first 
introduced in Ref. [42]. For our purposes it is useful to adopt a non-perturbative generalization of ( given by (see 
Refs. [20, 35, 43]) 

^ = -*+3(rT^i-?' ^'-''^ 

where just for simplicity in this expression we have assumed a constant equation of state parameter w = p/j5, p and 
p being the background energy density and pressure respectively. This variable lias two crucial features: it is gauge 
invariant and, most important, it is constant in time for adiabatic, super-horizon perturbations.^ This allows us to 
set the initial conditions at the time when ( becomes constant, and follow them until the perturbation mode reenters 
the horizon. The primordial non-Gaussianity can be parametrized in terms of the curvature perturbation as [20] 

C(^) = 2(aNL - 1) (C^^))' , C^'^ = 6(6nl - 1) {C^'^y , (4.18) 

where the two non-linearity parameters depend on the physics of a given scenario of generation of the perturbations. 
For example for standard single-field models of inflation onl = 1 and 6nl = 1 (plus tiny corrections proportional to 
the slow-roll parameters), while for other scenarios they might well be non- negligible. 

On the other hand, the physical observable quantity is given by the CMB anisotropics, which are one of the best 
tool to detect or constrain the primordial non-Gaussianity generated on large scales [3]. The standard procedure 
is to introduce the quadratic and cubic non-linearity parameters, /nl and gNL (which can be in fact also non- 
trivial kernels in Fourier and harmonic space) characterizing the non-Gaussianity in the large scale temperature 
anisotropics [1, 10, 14, 18, 44]. In the limit of large non-Gaussianity, |onl| '> 1 and |6nl| ^ 1 [3, 20, 45], the 
resulting size of the non-Gaussianity in the CMB anisotropics can be estimated as /nl ^ 5aNL/3 and ^nl ^ 25&nl/9 
if accounting only for the contribution to the CMB anisotropies from the Sachs- Wolfe effect. In fact for the cubic 
non-linearities entering in the Sachs-Wolfe effect such an estimate is too rough, and the correct expression relating 
the observable quantity ^nl to the primordial non-linearity parameters Cnl and 6nl is given by Eq. (67) of Ref. [20] . 

Our results allow to take into account also the contribution from the (late) Integrated Sachs- Wolfe eff'ect. The 
memory of the initial non-Gaussianity is kept in the matter-dominated value of the large-scale gravitational potential 
ipmK and also, as noticed before, in the source term (4.2). In order to determine iprn we use the expression (4.17) for 
the curvature perturbation in the matter dominated era 

Cm = -*m-^$m, (4.19) 

where we have employed the energy constraint (Bl) to replace for the energy density. 



^ As a consistency test, we have used Eq. (4.17) together with the Einstein equations contained in the Appendix, to verify that C' turns 
out to be proportional to the L.H.S. of Eq. (4.1) minus the very same source term appearing in Eq. (4.2). This provides an independent 
check for the expression of the source term which turns out to be consistent with ^' = on large scales. 
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We can expand both sides: at first order we find the usual result 



(4.20) 



at second order 



and at third order 



-3(aNL-l)Vo + Vo + ao 



With the use of Eq. (4.14) evaluated in the matter dominated period, we get as initial conditions 



r TO "Wto 



/^25„ 10, 2\ , „ 20, 

( y (^NL - 1) - y (oiNL - 1) + 3 1 "^0 + 2<PoQ!o - y (oinl - 1)Mo + 

25,, ^ 10, , 2\ , 4 160, 

y (onl - 1) + y (ONL - 1) + ^ I lySo " g^'oOo ^(onl - 1)Mo + ^vq 



Finally, from Eqs. (4.12), (4.14), the large-scale solutions for tp^^\ (f>^^^ are given by 



(4.21) 
(4.22) 

(4.23) 

(4.24) 



25„ 10, ,x 2\ o, , 4 , , , , 160, , , „ , ^ 

y (&NL - 1) + y (fflNL - 1) + - 1 </?o W - -<po(x)ao(x) - — (onl - 1)mo(x) + 8uo{x) 



T^o [fo+ m 



5(7?) / dfia{fi)n{fi)S{fi,^) . 



dr/a {fi)g{fi)S{fi,yi) 



0(-^)(r,,x) = V(')(r?,x) + 12.9(77) A(7?)^3 ^g^(^^ (^(77) + C(77)) ao(x)^o(x) 



(4.25) 
(4.26) 



V. CONCLUSIONS 

In this paper we have computed the expression for the CMB anisotropies due to the redshift the photons suiTer 
when they travel from the last scattering surface to the observer up to third-order in the gravitational perturbations. 
We thus have completed the findings of Ref. [20] , where a fully non-linear expresssion for the Sachs- Wolfe effect has 
been obtained, by including also the Integrated Sachs- Wolfe effect and lensing effects at the last scattering surface. To 
achieve this goal we have proposed an alternative method to the standard perturbative one in order to solve for the 
geodesic photon equation which provides a fully non-linear integral solution. Moreover we have studied the evolution of 
the gravitational potentials on large scales, allowing for generic non-Gaussian initial conditions. Our results, together 
with those of Ref. [20], are of particularly relevance when facing the trispectrum of the CMB anisotropies, in that, 
as pointed out in the Introduction, they include all the relevant cubic non-linearities for a coherent prediction of 
such statistic from various cosmological scenarios. In particular they include those contributions related to primordial 
non-Gaussianity, some of which would be inevitably lost if one would stick to the linear evolution approximation that 
is often used in the literature. 
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APPENDIX A: CONNECTION COEFFICIENTS 



We start from the line element of the conformal metric 



^drf + 2ujidr]dx' + (e"^*(5y + Xij) dx'dx^ 



It is useful to consider as "background" metric the one with the exponentiated scalars: 

900 = -e , goi=0, Qij = e dij , 
while treating vectors and tensors perturbatively (so here n refers to the order in vectors and tensors): 



(n) 
00 



0, 



1 



in) 



^9if = Ax 



in) 



n! ' ' nr'^ 
The inverse metric will be computed perturbatively up to third order. At zeroth order 



-00 



-2* 



5°^ = 0, 



gij = g2*jij 



At first order, Sg^^-. = —g'^^g^'^Sg'^^ and we find 



9», = p2(*-$)^i 



"5(1) - e X(i; 



At second order, Sg^^^ = -g^'^g^PSg^^], - Sg^^^-^g^PSg^^l and we find 



6gl% =e^*-^* (4^y 



At third order, Sg'^^^ = -g'^^g'^^Sg^^^ - (5; 



, ,i p2* (1) fei 

'^(2) - e Wj. X(i) 



,00 _^2*-4* 



'(1) 



^Ap- 



es) 



i ,.,(2) _ .2* 



„2(*-*) 



(3) 



5'^(3) ~ 2 * ^(2)~^ M,^; ) '^(i)~X(i)e 
6^(3) "2 (,'^(i)'^(2) + '^(2)'^(i) J + 2'' iX(i)X; +X(2)Xi ) 

„2(*-*), ,(1) / /i , J I , ,J A p4* ife ^(1) H 



e X(i) 



It is now convenient to resum the inverse metric, retaining terms up to third order; we find 

500^_^-2* + g2*-4* (^,)2_^Vxi,-^ 

gOi =e2(*-*)a;' - e^*-2*c.,x'^ " M'u;^ + x''x,kOJ^ , 

gij =e2*^i,- _ ^4^^ij _ e4*-2*^i^.- + e'^^x'^xl + {u^'x'' + ^'x"') - x'^Xkix'' ■ 
It is now immediate to compute the connections coefficients with the the usual formula 



Kf} = of {dagi3\ + d^gax - Oxgan) , 



where g^^"^ is the resummed inverse metric in eq.(A8). We find, keeping terms up to third order: 



r°- =\E''(5- 



+ 2 i^iXjk + djXik - dkXij) 



2*, k m 



rj,o = [-6^*0;* + ujjx''] + + e'*a,-$) [6^*^*^ - e^*x'^' - w'o;^' + x^'^k 



n^.=a,$[-62*o;'+o;fex'1 

+ 1 (x^.fe - 2*'6-2*5,., + a,o;fc - dkWj) [62*5'^= - e4*x'' - ^'^^ + x''xf] ; 



+ 2 [(^iXfez + dkXji - dixok) - 2e-2* {d,-^5ki + 9^*^,, - 
X [e^^S"' - e4*x'' - e**-2*a;^o;' + e6*x'"x^] • 



APPENDIX B: EINSTEIN EQUATIONS FOR THE GRAVITATIONAL POTENTIALS 

In this Appendix we derive the evolution equations for the third order scalar perturbations (l>^^\ 
We start by writing the fully non-linear Einstein equations for the metric ds"^ = a'^{r}) \—e'^'^drf + e-^*^,^ 
for a ACDM model: 

• 0-0 

6-2* (T^-M/f = ^(87rGp + A) ; 



• 0-i 



6-* [di^' + {n- = -AwGapui ; 



• i-j traceless 



d'dj^ - -V^^S'j - d'dj^ + -V2$ -I- d'^dj^ - -d^^di'^S'j 



=87rGa26-2* ( Tj - ^5}T 



• i-j trace 



-2* 



[{H - {-H + 3*' + 2$') - 2H' + 2*"] 



„2* 



[2V2 ($-*)-!- -h 2a'$5i$ - 25'$5i\E'] = 



SttG 2 
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In writing down the first two equations, we made use of the normalization condition for the velocity, g^^u^Ui, = —1, 
which gives 

uo = -e* (a^ + e2*u'«i) ' ~ -ae* ; (B5) 

the last equality is valid at any order on large scales, and up to third order on smaller scales. 

We start from the traceless equation, which we project onto the scalar modes by applying the operator did^: 



(B6) 



It is convenient to rewrite it as 

* - $ = Q , (B7) 

where we define 

V'^Q = -P + 3N = -PI + SV-'^did^P] (B8) 

and 

P; = 1 [{d'^M^) - {d'mdj^ + \ + + 47rGa2e-2*r] . (B9) 

For a ACDM model, the spatial part of the stress-energy tensor is 

. e2* - A - 



a2 ^ SttG ^ 

^ *^ [d'^' + {H- [dj^' + {n- dj^ - -^s] 



(BIO) 



167r2G2 aV ^ ^ ' ' SttG ^ ' 

where in the last equality we substitute for the spatial velocities using the — i equation: 

Ui = --^—m' + {n-^')di^. (Bii) 

Substituting this expression into we find 

[5**' + {n- d'^] [a,*' + {n- a^^] - -e-2*a2A(5i . 



(B12) 



47rGa2 p l ' v / j l ^ > v ^ ^ 2 

Take now the trace equation, and using the background equation + 2H' = a^A and eqs. (B7), (B8), we find 

[-a^A - 2*'$' - 3(*')^ + 2W(3*' - Q') + 2*"] 

+ — [2d^<^di^ + + -25'$5;* + 2(P - 3N)] = -—a^T . 

o o 

At this point, we can expand eq. (B13) at third order. Using the first order solution = "^^^^ = cp, we find: 

^(3)" + 37^^(3)' + a^A^f^) = 6a2AQ(3) + Sa^Acp^f^) _ ^^2^^3 _ 3o^(<^')2 + 3^' (4^(2)' + $(2)') 

+ ISW^^^V - 36WvV + 18W(^*(2^' - 12W(^Q(2)' ^ eWQ^^^' 

where Q = Q'^) + Q(2) + qO) _^ _ _ _ _ 
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Since wc are interested in the large scale solution, it is more convenient to use the small letter variables tp^'^\ 4'^^\ 
related to the capital letter ones as explained after eq. (2.6). Finally, using the first order equation (p" -{-'iTi^p' +0^ h.^p = 
0, we get 



At this point we need to explicit the terms Q^^\ Q^^^' and T^^\ It is straightforward to show that 



(B15) 



J.(3) ^ 



(47rGa2)2po 



and 



(B16) 



(B17) 



- 2np' d'pdjp + d^p' + av' a^-^f^) + d'^^^'djp + d'pdj¥^'^'j 

- p>' [5V' djV + d'pdip'\ - ((5(1) + 2(^) (aV' + ipi'p' + M'^') } • 

We now use the solutions at first and second order, and the first order (0-0) equation (5^1) = —2-S^ipo, to obtain 

+ F(r?) [2,V-^did^ {d'ao dj^o + d'p^odjao) - 2{d''ao dk^o)] } , 

where E{r]), F{ri), ao(x) are defined in eqs. (4.10), (4.11), (4.6) respectively. 
By differentiating the last equation we find also 



(B18) 



1 



{E' [3V-^did^ (^0 aVo dj^o) - (^0 aVo dkipo)] 



F' [3V-^did^ (aVo d,po + d'podjao) - 2{d''ao 9fe^o)] } 



(B19) 



using the background equation pQ = —3Hpo; and 



y(3) ^ 



ig' + 'Hg) 2A' + nA + nC- 2gg' + 2n 



(47rGa2)2po 
+ {g' + ng)B' {d'^aodkPo 



g^ 



(B20) 



where A{r]), B{ri), C{rj) arc defined in eqs. (4.7), (4.8), (4.9) respectively. 

Finally, putting all these expressions together, we arrive at eq. (4.1), while eq. (4.14) is simply given by expanding 
eq. (B7). 
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